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1 Introduction 

In this paper we study the mean-field Potts model with Hopfield-Mattis 
disorder, in particular with Gaussianly distributed disorder. This model is a 
generalization of the model studied in ||BvEN | . It provides yet another exam- 
ple of a disordered model with infinitely many low-temperature pure states, 
such as is sometimes believed to be typical for spin-glasses [[MP V | . In our 
model, however, in contrast to [BvEN], instead of chaotic pairs we find that 
the chaotic size dependence is realized by chaotic q(q — l)-tuples. For the 
notion of chaotic size dependence, and the notion of chaotic pairs which were 
introduced by Newman and Stein we refer to |N3L [NS2|| and references men- 



tioned there. Compare also [BvEN| and ||Nid] . For an extensive discussion of 



the Hopfield model, including some history and its relation with the theory 
of neural networks, see [[B], pag. 133 and further] or ||BG|| . A somewhat dif- 
ferent generalization of the Hopfield model to Potts spins can be found in . 

We are concerned in particular with the infinite-volume limit behaviour of 
the Gibbs and ground state measures. The possible limit points are labeled 
as the minima of an appropriate mean-field (free) energy functional. These 
minima can be obtained as solutions of a suitable mean-field equation. These 
minima lie on the minimal-free-energy surface, which is a m(q — l)-sphere in 
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the (e 1 , • • • , e q )® m space. This space for g-state Potts spins and m patterns 
is formed by the m-fold product of the hyperplane spanned by the end points 
of the unit vectors e q , which are the possible values of the spins. But only 
a limited area of the minimal free energy surface is accessible. Only those 
values for which certain mean-field equations hold, are allowed. These equa- 
tions have the structure of fixed point equations. We derive them in chapter 
H To obtain the Gibbs states we need to find the solutions of these equations 
on the minimal free energy surface. 

The structure of the ground or Gibbs states for q = 2 and £ fc Gaussian with 



m = 2 is known since a few years [ByEN]. Due to the Gaussian distribution 



we have a nice symmetric structure: the ground states form a circle. For 
a fixed configuration and a large finite volume the possible order-parameter 
values become close to two diametrical points (which ones depend on the 
volume of the system) on this circle. This paper treats the generalization of 
this structure to g-state Potts spins with q > 2. To have a concrete example, 
we concentrate on the case q = 3. It turns out that we again obtain a circle 
symmetry but also a discrete symmetry, which generalizes the one for Ising 
spins. One gets instead of a single pair a triple of pairs (living on 3 separate 
circles), where for each pair one has a similar structure as for the single pair 
for q = 2. For q > 3 we get g( - g ~ 1 - > pairs and a similar higher-dimensional 
structure. 



Our model displays quenched disorder. This means that we look at a fixed, 
particular realization of the patterns. It turns out that there is some kind of 
self- averaging. The thermodynamic behaviour of the Hamiltonian is the same 
for almost every realization. This is the case for the free energy and the as- 
sociated fixed point equations, as is familiar from many quenched disordered 
models. However, this is not precisely true for the order parameters. We 
will see that they show a form of chaotic size dependence, i.e. the behaviour 
strongly depends both on the chosen configuration and on the way one takes 
the infinite- volume limit iV —>■ oo (that is, along which subsequence). 



2 Notations and definitions 

We start with some definitions. Consider the set A;v = {1, ■ --,N} G N + . 
Let the single-spin space % be a finite set and the N-spin configuration space 
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be x® N ■ We denote a spin configuration by a and its value at site i by t7j. 
We will consider Potts spins, in the Wu representation |[Wu|1 . The set is 
then the iV-fold tensor product of the set \ = { el ? ' ' ' > e<1 }- The e CT are the 
projection of the spinvectors e a on the hypertetrahedron in R 9_1 spanned 
by the end points of e a . For q = 3 we get for example for e 1 , e 2 and e 3 the 
vectors: 

f / 1 \ / 
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The Hamiltonian of our model is defined as follows: 







m N 



-P H * = §EE &$Kws\ with 



k=l i,j=l 



where is the z-th component of the random iV-component vector £ k . For 
the £f's we choose i.i.d. N(0, 1) distributions. The vectors £, k = (£*, • ■ ■ , 
by analogy with the standard Hopfield model, are called patterns. If we 
combine the above, we can rewrite the Hamiltonian as: 



-I3H N = (3 



So asymptotically 



k=l 



Egg 
N 



+ 
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with K = 2(3 



g-1 

, 9 . 



iV 



and order parameters g^jv = T7 X! £f e<J ' : 



iV 



8=1 



(1) 



The last term is an irrelevant constant; in fact it approaches zero, due to 
the strong law of large numbers. (The £f's are i.i.d. N(0, 1) distributed so 
= 0.) Note that any i.i.d. distribution with zero mean, finite variance 
and symmetrically distributed around zero will give an analogous form of 
H N , but we plan to consider only Gaussian distributions, for which we find 



that a continous symmetry can be stochastically broken, just as in [[BvEN 



From now on we drop the subscript N to simplify the notation, when no 
confusion can arise. 
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Furthermore we introduce two representations for the order parameters q. 
If we assume m = 2 then q = (q\, q 2 ) and the definitions are as follows: if we 
consider the space R ql spanned by the vectors e 1 , ■ • ■ , e q , the x-plane, we 
define q = (x\, ■ ■ ■ , x 2 ( q -i))- It is often more convenient to look at the (higher- 
dimensional) (e 1; • ■ ■ , e g )-space. In that case we take q = (a 1; a 2 , a 3 , 6 1; b 2 , b 3 ) 
for q = 3 and an equivalent equation for other values of q. For m 7^ 2 the 
definitions are analogous. 



3 Ground states 

Now it is time to reveal the characteristics of the ground states for the Potts 
model. First we discuss the simple behaviour for 1 pattern. Then the more 
interesting part: q > 2 and 2 patterns. 



3.1 Ground states for 1 pattern 

For one pattern £ the Hamiltonian is of the following form: 

K B N 

We easily see that the ground states are obtained by directing the spins with 
£i > in one direction and the spins with ^ < in a different direction. 
If we have as the distribution for the £$'s P(£j = ±1) = \ , then the order 
parameter is of the form: q\ = |(e°" 4 — e"^ ), with 1 < i, j < q and i ^ j, 



see also ||vEHP|| . So for q = 3 we have only 6 ground states. They form a 
regular hexagon: (±|,T^)j 4); (0|i 2 )■ This regular hexagon with 

its interior is the convex set of possible order parameter values. It is easy to 
see that for N(0, redistributed we get the same ground states except for 



a scaling factor y2/n multiplying the values of the order parameter values. 

3.2 Ground states for 2 patterns 

The Hamiltonian for 2 patterns (Gaussian i.i.d.) is: 

3 N K 
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Similarly as in [|B vEJN | , we make use of the fact that the distribution of 2 



independent identically distributed Gaussians has a continous rotation sym- 
metry. This symmetry shows also up in the order parameters. Let 

fffa\-( x ^°) \-( as[n9 \ ff(o\-( x a( ) \ -( acos9 \ (9\ 

with (a, (3) a ground state associated to the special case 9 = i.e. to the 
second pattern. We note that asymptotically for large N we get the same 
ground state energy per site for each value of 9. Because the surface on which 
the Hamiltonian is constant is of the form gf + gf = C 2 , these are the only 
ground states. For finite N, however, there are finitely many (q(q—l)) ground 
states, corresponding to one particular value of 9 (The exact symmetry of 
choosing a different pair of Potts directions gives the g(g — 1) ground states). 
This is an example of chaotic size dependence, based on the breaking of a 
stochastic symmetry, of the same nature as in [ BvEN|| . Because of weak 
compactness, different subsequences exist whose q(q — l)-tuples of ground 
states converge to q(q — l)-tuples, associated to particular ^-values. These 
subsequences depend on the random pattern realization. See Appendix |A]. 
For further background on chaotic size dependence and its role in the theory 
of metastates we refer to |N S| | . For m > 3 patterns one has the same discrete 



structure as before, but instead of a continous circle symmetry we have a 
continous m-sphere symmetry (isomorpic to 0(m)). 



4 Positive temperatures 

At positive temperatures instead of minimizing an energy one needs to min- 
imize a free energy expression. 

By making use of arguments from large deviation theory we obtain (see e.g. 
-Pf{P) = sup{Q(?i, &) ~ c*{q!, q 2 )}, 

91,9*2 

where / is the free energy per spin and —(3H = A^y(gf + gf) = NQ. The 
function c* is the Legendre transform of c, where c is defined as follows: 

c(i) = lim — In (E CT exp(t! • iVgl + t 2 ■ Nq 2 )\. 
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Here t\ and t 2 are vectors in R 9_1 and tr CT is the normalised trace at a single 
site. To determine the supremum (maximum) we differentiate and put the 
derivative equal to 0. This implies that for q\ and q 2 it holds: 



?2)max = Vc(VQ(gi, <f 2 )) = Vc(Kq u Kq 2 ), with 



(3) 



Kq 2 



ag_ 

dm 
dQ 

9q2 ' 



We make use of the fact that for a convex function c, Vc* = ( Vc) 1 (see also 
|[BG| , chap. 3] and compare pag. 27). Now let us rewrite c(i): 

c(t) = lim — In (E CT exp(fi • Nq x + t 2 ■ Nq 2 )\ = 

• • • =< lntr CT {exp (^ + rjt 2 ) ■ e a } >^ v . 

Plugging this into @ we get the mean field equations for the order parame- 
ters which have the structure of a system of fixed point equations q = F(q): 



tr CT {ge CT cxp [K(ggi+i7g 2 )-e' T ]} 
tr CT {cxp [if(£gi+W2)-e <T ]} / £ >r] 
/ tr ; {? ? e' r {cxp [K(ijq 1 + m2 )-e' J Ys 
\ tr (7 {exp[A:(^gi+7 ?9 2)-e CT ]} 



(4) 



If we are in the allowed area, that is, the domain of definition of F, it is 
equivalent to look in the (ei, • • • , e q ) -space. We may rewrite @ in this area 
as follows: 



Qi 



\ 



Q2 



\a q / 
(h \ 



I 



gexpJC(^ai+7?bi) 



\ 



£ cxp K (£a q +r]b q ) 



( 



r) exp K^ai+rjbi ) 
J2i=i exp K (£ a i+Vbi) / ^ 

rj cxp K (tjctg+ribq ) 



\ 



\ \i:Li ex p^(^+*)/ e , r) / 



with gi = E?=i a^i and q 2 = Ya=i h&i- 
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4.1 Ising spins 

If we look at the behaviour for iV —>■ oo, then due to the strong law of large 
numbers 4Si=i6 = E£ = 0. Each coordinate aj of vector q\ = (ai,a 2 ) 
is defined as jr Y^Li Ci^i^h This means that aj is the contribution of 
the spins in the j-th direction to the sum jrJ2iLi£i- Therefore: a\ + a 2 = 
jr J2i=i & = a - e - This gives a necessary condition for the allowed area of 
Ising spins: 

ai = —02 A b\ = —b 2 . (5) 
Furthermore for all Gibbs states the value of the energy is constant, therefore: 

a\ + al + b\ + bl = — . (6) 

When we substitute in equation @ and project the result to the (xi,x 2 )- 
plane by the projection II : ei — > 1, e 2 — > — 1, we obtain the following 
equation: 

2,2 *2 

x 1 + x 2 = r . 

Thus to get the radius of the circle of Gibbs states r*, just take the point 
<fi = (a, —a), <f 2 = (0, 0). This corresponds to the point (2a, 0) in the (xi, x 2 )- 
plane, by the projection IT. Of course 2a = r*. 

With this we calculate the equation for the first coordinate of g! in the 
(ei, e 2 )-plane by substituting the corresponding fixed point equation: 

2tt 7 7 exp/3£a + exp(-/3£a) \ ° 1 



exp /3£a 

exp 




\[2tx J exp /3£a + exp (— f3£a) 

We replaced .K" by (3, because for Ising spins K = 2/3(2 — l)/2 = (3. The 
equation for the second coordinate of q\ we calculate in the same way. The 
vector q 2 is simply (0,0). Now project q\ and <f 2 to the (xi, a; 2 )-plane. That 
is done by subtracting the second coordinate of the q^s from the first one. 
We get the following equation for the radius r*: 

r* = -L / etanh ( ^] exp f-^W (7) 



'2ti 
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For (3 > /3q this equation has a nontrivial solution for r*. The equation is 
the same as in jByEN] except the factor 1/2 in the tanh. This is due to our 
using the Wu representation. 



4.2 Potts spins 

If we take q = 3, then K = The set of ground states now can be 



parametrized by three (in general 



g(g-i)' 



circles, and similarly for the low- 



temperature Gibbs states. To obtain the radius f of such a circle parametriz- 
ing the ground or Gibbs states, we follow the same recipe as in the case of 
Ising spins. Here we take the point (<fi,<f 2 ) with q\ = (0,r/y/3, — r/v3) and 
<?2 = (0,0,0) (the representatives of both qi in the (ei, e 2 , e 3 )-plane). Now 
(fi projects to (0, f) by the projection: 




_ i 
£V3 




So if we substitute the corresponding fixed point equations for q\ in the 
(e 1; e 2 , e 3 )-plane, we get for the order parameter values (ai, a 2 , a 3 ) = q\ the 
following mean field equations: 






I r £ exp (K 

75^ ' ^ exp (K£r/V3)+cxp (-K£r/\/3)+l 

exp (-Kt-f/y/3) 

exp(A'gr/v%)+exp(-X5f/v / 3)+l 
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Here (a 1; a 2 , a 3 ) = (0, r/v3, — r/v3). Thus by taking the difference between 
a 2 and a 3 and multiplying it by |a/3, we finally get the following expression 
for the absolute value f: 




exp (K£r/V3) ~ exp {-K^f/y/^) 
exp (K£r/V$) + exp (-K£r/V3) + 



exp 





£sinh (K£r/V$) 




2 cosh {K£r/V3) + 1 6Xp 1 ° 1 d4 ' (8) 
We can easily check that this expression indeed approaches the one for the 
radius for the circles through the ground states, by considering the behaviour 
of the integrand for K — > oo. It behaves like: 



/i«i«p(4)* 
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A Stochastic symmetry breaking for q = 3 



In this Appendix we adapt the fluctuation analysis of PvEN|| to include 
Potts spins. We essentially follow the same line of argument, and find that 
the fluctuations, properly scaled, after dividing out the discrete symmetry, 
approach again a Gaussian process on the circle. 

For notational simplicity we treat the case q = 3 only. For q > 3 a simi- 
lar analysis applies. Define the function <pN,2 as follows: 

P(t>N,2{z) = -Q(z) + z- VQ(z) - c(VQ(z)), 

where c(i) equals: 

1 1 N 

c(t) = — In {E CT exp h ■ Nq\ + t 2 - Nq 2 ) = — ^ In {e ct , exp h ■ &e CT * + t 2 ■ Vi e^}. 

This is chosen such that for A" — > oo the measure 

p -f3N<t>N 

C= e — >£, 

where £ is the induced distribution of the overlap parameters. 
For q = 3 it holds: 

QW = ^\\m=bm\ 2 - 



Thus: 

2 1 r 4 1 2 1 _ 

V.2O - In {E a exp ■ e^+r)^ ■ e ai ) \ = -\\z\\l-—E N<2 . 



3" " z f3N I " r 3^ v * 1 ' ') 3" iU (3N 

_ 2 K . KVS,, 

-n,2 = 2_, ln j 3 ex P K \&\i + ViZn) + 3 ex P -y (6^11 + ^21) cosh — ^ — (6^12 + ^22) 



N 



1 , , . 2 



H ln < 1 g^i(^ii) «22){,i, + zrn=i ^=M z i2,z 22 )^ v 



i=l 



3^J 0i(zn, z 22 )^ 



Because for finite A" the set of 6 Gibbs states has a discrete symmetry, as 
mentioned before, we choose out of these 6 states one state we like, namely 
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the one of the form (0, iasin^, 0, ±acos6*). Note that the 9 depends both 
on N and on the realization of the random disorder variable. Then z n = 
-^21 = 0i = 0. Inserting this and defining Z\ 2 = Z\ and z 22 = z 2 we get for 0: 



Putting (zi,z 2 ) = 3$(zi,z 2 ) we obtain: 



1 f 1 2 2 

zi, z 2 ) = - 1| (ii, z 2 )\\l - Yl ln 1 3 + 3 cosh ^ P(£& + ^2) 



1 1 N ( 1 112 

0Oi,2 2 ) = 2 il(^i^2)||2-^ I] In |- + cosh^(^i+^ 2 )| ~^] v ;ln 3- 

^From now on the last term will be ignored. So it is enough to prove now 
that with the | term we get the desired chaotic pairs structure between the 
patterns due to the quenched disorder for this class of ground states, once 
we divide out the appropriate discrete Potts permutation symmetry. Thus 
the original model displays chaotic 6-tuples. 

Therefore we only need to control the fluctuations of 0. Define 

I N 1 N 

r N (z)-Er N (z) = — 5>{cosh0?. (e,^)}- — ^Eln{cosh^- (£,77)}. 

(9) 

This is the fluctuation of the Ising case which we can estimate by |[BvEN|| . 
Denote the corresponding function by 0*. We start with the following 
lemma: 

Lemma A.l 

exp (-(3N(j)) < exp (-/3iV0*). (10) 

Proof: 
Because 

exp (-(3N(j)) = exp (-/3iVE0*) exp {-/3N(<f> - E0*), 

we only have to estimate the quantity 0— E0*. Notice that also a lower bound 
is essential, because the quantity can become negative. First the estimate 
from above: 

- E0* = — £ In I - + cosh ■ (£, V ) I - — £ E ln {cosh (3z ■ (t, v ) }. 
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Now use 



In { \ + C08 h Bz ■ ({, „)} = In {l + ^^(^ } + 'n {«»> ,3? .({,,)}< 



In {cosh f3z ■ (£, 77)} + In - to get 



1 N 1 iV IS 

)— E0* < — glBfoosh^. (6,77)}-— EEln{cosh/3z- (£,,,)}+- In- 



A? 



y*(3-E&(3 + ±ln| 



(11) 



This is because cosh x > 1 for all x G R. 



The lower bound is easy because: 

1 N rl i l * 

This is due to the fact that the function In a is monotonically increasing in 
a. Then it follows that: 



^>r N {z)-^r N {z). 



(12) 



Combine (JTTJ) , ([12]) and use the fact that in the limit lim^^co the constant 
term jj In | does not contribute to the expression exp {— pW(0 — E0*)} to 
conclude the proof of lemma |A.1| . 



Henceforth it is convenient to transform <ft* to polar coordinates. Define 
z = (r cos 8, r sin 6). Then (O) transforms to: 



\fUr,0)\ 



N 



1 1 

-E^E, : In cosh {(3(r cos ^} - -jtj m cosn {^ r C* cos _ ^i) } 



|E&(r,0)-/&(r,0)|. 

Here £, V* denote the polar decomposition of the two-dimensional vector (£, 77), 
i.e. ( is distributed with density xexp —x 2 /2 on R + and ip uniformly on the 
circle [0,27r). See [ |BvKl\l| , page 188]. This we see easily because: 



£zi+rjZ2 = (( cos i/j) (r cos 6) + (( sin i[))(r sin 6) = £r(cos#cos'?/>-|-sin6 l sin'?/>) 
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(r cos {8 — i/j) and E^, cos {8 — ip) = E^, cos if). 
With 0* in this form, estimate (|TDD of lemma [A.l| is not very useful, since 



the fluctuations of reach their minimum for a different radius (in f) in 
general than the fluctuations of 0* (in r*). Thus we need to transform 0* 
such that the fluctuations of the transformed 0* reach their minimum at the 
same radius f as those of 0. This we achieve as follows. There is a uniform 
transformation II which translates all the points on the circle with radius r* 
centered at the origin to the circle centered at the origin with radius f, the 
radius of 0. If we apply II to 0*(r, 8) then we get 0*(r + r* — f), the desired 
transformation of 0*(r, 8). Now we can prove the next lemma: 

Lemma A. 2 For every e > holds: 

\fN(r,8)\ < |/*( r + r *-f,0)|+e. (13) 

The constant r* is the radius of the circle parametrizing the set of mean-field 
solutions in the Ising case (q = 2). The constant f is the radius f in the 
Potts case q = 3 rescaled by the factor 2/ a/3, thus f = (2/ y/3)r. 

Proof: 



We use the following estimate, which is lemma 2.5 from [BvEN 



|E/^(r, 8) - fx{r,9)\ < | a.e. on every bounded set. (14) 

Define: 

= {zeR 2 : \\z\\ >r* + 5}, a = {zeK 2 : \\z\\ >r + S}. 

Set O C O' because f < r* '. Check this by using (|?p and (^) and the scaling- 
factor 2/y/3 for f. Decompose O' as £> U C \ O. Because O' \ O is a finite 
set we can use estimate ([HI). With the already obtained estimate for O in 
BvEN] , (|13D holds for all (r, 8) G C. Because of (0) is true for all finite 



sets, (0) holds for all (r,6). 

Note that in a neighbourhood of f it is equivalent to look in a neighbourhood 
of r*. Now we are able to prove the following theorem: 
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Theorem A.l Let £ be the induced distribution of the overlap parameters 
and let m = m(8) = (r cos 8, r sin 8), where 8 e [0, 7r) is a uniformly dis- 
tributed random variable. Then: 

Furthermore, the (induced) AW-metastate is the image of the uniform distri- 
bution of 8 under the measure-valued map 8 — > C 00j p\m{8)\. 

First we prove the following two lemmas: 

Lemma A. 3 For <f)^ one? 7] i; with i e N as defined above, there exist 
strictly positive constants W, W, I, I' such that (f is the largest solution of 

®) 

J\\\z\\-f\>5 N e of Wp~ wn1 

i i 1 i 
on a set of V -measure at least 1 — W'e~ K N , where 5n = N~w. 



Lemma A. 4 Assume the hypotheses of lemma A.S . Let a N = 



len 

there exist strictly positive constants K\,Ki,C\,Ci such that on a set ofP- 
measure at least 1 — K\e~ Nl! ^ the following bound holds, 

j -pN4>„w d z Arl/5 

JA N < p p-N 1 ' 5 

J /if] 

where 

A N = {(r,8) eR + x [0, 2tt)| |r - r\ < 8 N , g N (8) - min e g N (8) < a N } 
A' N = {(r,8) G Rq x [0,27r)||r-f| < 5 N , g N {8) - min e g N {8) > a N }. 
Here 



9n{8) 



n r l i l r l 

E^,E^ In < - + cosh (3(r cosi/j > -= 2J m l o + cosa {P r d cos 1 



P 

which is the polar coordinate form of the function gN{z), which is defined as: 
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It is convenient to look at the following decomposition: 

(0jv - E<j) N )(z) = (3VN{g N {2) + h N (z)), where 

h N (z) = g N {z) - g N {z). 

The variable z 1 is the projection of z onto S' 1 (f). Note that (3y/Ng^ = 
4>n — E0iv = In- Define g* N and h* N in the same way but as decomposition 
of instead of Jn- 



Proof of lemma \A.3{: 



Compare lemma 2.1 in |BvEN| 1 . Define: 

= {zeR 2 : \\z\\ >r* + 5}, a = {zeR 2 : \\z\\ >r + S}, 

1 = {zeR 2 : \\z\\ <r*-5}, I' = {z GR 2 : \\z\\ <r-6}. 
Now we first estimate the numerator which we can also write as: 



[„ 

■/III 



z\\— f |>(5jv 



-/3N<f> N (z) d g 



e -(3NV<j> N (z) e -f3N{<i>N (O-E^iv (i*)) _ 



O'VJV 



With lemma [A.2| we have the following inequality: 



sup |/&(r,0)| - e < sup \ f* N (r + r* - f, 0)\ = sup |/^(r) 
zee zeo' zao 



sup 

,(r,6»)eC 



c 

e > — (r — f) 2 
- 2 v ) 



< P 



sup \r N ( r ,6)\>^(r-r*) 2 . 

(15) 

Lemma 2.4 of |BvEN|l tells us that this event is of measure zero. Now we can 
estimate the integral. 

First we estimate E0^(,?). Because F,(p* N (z) is a bounded function, in each 
bounded interval one can always bound it from below by a function of the 
following kind: 

~E(fi* N (z) > C( 1 1 z\ | — f) 2 + E0^(f), with C a positive bounded constant. 
Then 

E(j)* N (z + r* - r) > C(\\z\\ - r*) 2 + E^(r*), 
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when we apply II to this estimate. Now use estimate ( |To] ) with this constant 
C. Then it holds: 



/ e -/3irE4,irW e -l3N(4>irW-B4>N&)d z < / e -/3ArE^(?+r*-f) e /37V|/* 0OI e ^iV rf - < 
JO' ~ Jo ~ 

e -/3iV(E^(r*)-e) f e -f3NC(r-r*) 2 e f3N%{r-r*f ^ _ e /3JV(E^(r*)- £ ) /" g - /37V §(r-r*) 2 rfr < 

" - 27r ^ exp (-w^mo - e )) ex P -p nc * (j) ■ 



For further details see PvEN|| . The interior X gives a similar expression. 



Notice that the image of X under the transformation II: r — > r + r* — f is 
X \ 73(0, r* — f). The ball -8(0, r* — f) is a finite set so we can integrate 
over X instead of X \ -8(0, r* — r) by ([14]). To estimate the denominator we 
just replace r* by f in the expressions of the proof of lemma 2.1 in |[BvElN 



pag.192,193]. Combining the estimates for the numerator and the denomi- 
nator gives the desired result. 



Proof of lemma \A.4 



^From this moment we ignore the constants which enter by applying lemma 
A. 1|, because they cancel out when we divide the numerator by the denomi- 



nator. For |/ijvi it holds: 

l^/vl < \h* N \ < e, 



by lemma 2.6 of [BvEN ] . Consider the following integral: 



J8:g N (e)>a N +mm e g N (8) 



Ng N (e) < 2 7re -v / Jvmin e9JV (6») (r V v-,^ 



Henceforth it is just a matter of plugging in to get the desired estimate on 
the denominator. We refer to IBvEN] for the details. One gets a estimate 
for the denominator in the same way. By dividing the two estimates lemma 



is proven. 



Proof of theorem \A.l 



In the preceding paragraphs we have seen that the measures C concentrate 
on a circle at the places where the random function gjv(#) takes its minimum. 
Now it only remains to show that these sets degenerate to a single point, a.s. 
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in the limit iV — > oo. If we have proven it for L, then we have proven it also 



for C, because liniAr^oo^ = C. With the help of [ BvEN | this is very easy. 



because we can use Proposition 3.4 with the function 

g{.) =ln|cosh / S. + -|. 

This works because g is an aperiodic even function. And of course Proposition 
3.7 also holds for this g. These two propositions we use, tell us that the 
process r] N = g7v(6>) — ~Eg N (9) converges to a strictly stationary Gaussian 
process, having a.s. continuously different iable sample paths. And on any 
interval [s,s + t], t < tt the function t]n has only one global minimum. 
Furthermore, if we define the sets: 

L N = {9e [0,7r) : rj N (6) - min^^') < e^}, 
with €n some sequence converging to zero, — > 0*. Then the remarks 



below Proof of theorem 3 in [ BvEN ] conclude the proof. 
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